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Abstract -We present calculations of Majorana edge modes in cylindrical nanowires of a semicon- 
ductor material with proximity-induced superconductivity. We consider a Rashba field along the 
transverse direction and an applied magnetic field in arbitrary orientation. Our analysis is based 
on exact numerical diagonalizations for the finite cylinder and on the complex band structure for 
the semi-infinite one. Orbital effects are responsible for a strong anisotropy of the critical field 
for which the effective gap vanishes. Robust Majorana modes are induced by the parallel field 
component and we find regimes with more than one Majorana mode on the same edge. Exper- 
imentally, they would manifest as a specific sequence of zero-bias conductances as a function of 
magnetic field. In the finite cylinder, a degradation of the Majorana modes due to interference of 
the two edges leads to oscillating non zero energies for large enough fields. 



Introduction. — Novel states of matter, such as 
quantum Hall systems [HI], obey topological rules. The 
advent of topological materials [5HH] was a major break- 
through towards the achievement of a topological quan- 
tum computer. In these materials, quasiparticles are non 
Abelian composite fermions [9Tlllj. the potential basic 
units for quantum information processing [T2]- Such pe- 
culiar excitations are identical to their own antiparticles 
and they are termed Majorana fermions (MFs) honouring 
Ettore Majorana who predicted their existence in 1937. 

Majorana's original aim was to explain the nature of 
neutrinos as real solutions of the Dirac equation for rel- 
ativistic particles [13T[T7] . Notoriously, there is not yet 
a conclusive clue on the existence of MFs as elementary 
particles but they have been engineered in a laboratory as 
quasiparticles. Majorana fermions in topological materi- 
als can be described as half-fermionic states, so that a pair 
of MFs forms a widely distributed full-fermionic state. A 
pair of MFs defines a qubit whose information is encoded 
not in the individual particles but in a non local pair of 
identical and neutral entities. Non locality gives a great 
advantage to such pairs, they are immune against local 
sources of decohercncc. Majorana fermions are incredi- 
bly elusive because of their lack of characteristic features. 
They are chargeless, spinless and energy less, what makes 
them difficult to detect. 

The discovery of topological systems brought differ- 
ent proposals for the physical realization of MFs in solid 
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Fig. 1: Hollow semiconductor nanowire with cylindrical ge- 
ometry of length L. Cylindrical coordinates and unit vectors 
are indicated. A magnetic field B is applied in an arbitrary 
direction. 



state devices [T7H20] . Recently, they have been pro- 
posed in semiconductor nanowires as a combined effect 
of superconductivity, spin-orbit interaction and magnetic 
field j!?TM23| . Mourik and coworkers reported in Ref. {21] 
signatures of such quasiparticles in long InSb nanowires 
(L rj 2[im). In the experimental setup an InSb nanowire 
is connected to a normal-metal gold electrode on one side 
and to a superconducting niobium titanium nitride elec- 
trode on the other. Since the NbTiN contact is brought in 
close proximity to the InSb nanowire, leakage of Cooper 
pairs induces superconductivity, at least over a coherence 
length of a few hundred nanometers. InSb is character- 
ized by a large g factor, so that the application of an ex- 
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ternal magnetic field closes the superconductor gap, shift- 
ing some states to zero energy. The spin-orbit interaction 
induces anticrossings that separate the zero energy state 
from the nearest excited states, preventing the decoher- 
ence of the Majorana mode. The Majorana-mode evidence 
is a robust zero bias peak in the tunneling differential con- 
ductance, a zero-bias anomaly (ZBA), when the nanowire 
is magnetically driven into the topological phase. 

More experiments have reported evidence of Majorana 
modes in nanowires. Deng et al. have also observed a ZBA 
in the nonlinear transport measurements with an InSb 
nanowire quantum dot |25| . Almost simultaneously, Das et 
al. provided similar evidences in shorter hybrid semicon- 
ductor/superconductor InAs/Al nanowire junctions [26]. 
Rokhison and coworkers |27j have presented indirect mea- 
surements of the fractional a.c. Joscphson effect in hybrid 
InSb/Nb nanowire junctions as a hallmark of the topolog- 
ical matter. A more recent experiment by Finck et al. [28] 
demonstrates a striking effect: the splitting of the ZBA at 
high magnetic field and its oscillatory evolution, in agree- 
ment with the theoretical expectations [2TJ1 15D] , 

From the theoretical side, most works dealing with Ma- 
jorana physics in hybrid semiconductor/superconducting 
nanowires assume planar geometries. However, all the 
above mentioned experiments are done with cylindrical 
nanowires. This motivates us to investigate the theory 
of Majorana edge modes in cylinders. We quantify the 
importance of the cylindrical configuration, rather than 
a planar one, regarding the different influence of both 
Rashba field and orbital magnetic effects. 

We find that in cylinders a fixed Rashba direction has to 
be assumed, since a radial one docs not favour the forma- 
tion of Majoranas. Orbital effects are inevitably present, 
however, we demonstrate that they are less dramatic in 
cylindrical nanowires than in planar ones. In the latter, 
a slight deviation of the magnetic field from the nanowire 
plane completely destroys the Majorana state [55]. Re- 
markably, orbital effects drive the cylinder into new topo- 
logical phases in which several Majorana modes can co- 
exist on the same edge. We predict the specific sequence 
of topological phases in parallel field for finite and semi- 
infinite cylinders. Besides, we show that the critical field 
for gap closing strongly depends on the magnetic field ori- 
entation due to the presence of orbital effects. 

Model. — We study the spectrum of a semiconductor 
nanowire with cylindrical shape in a magnetic field. The 
quasiparticles are assumed to move only on the surface of 
the cylinder, with a total Hamiltonian H = Hkin + Hr + 
Hz+Hpair, where the successive contributions are kinetic, 
Rashba spin-orbit, Zeeman, and pairing ones. We include 
the orbital effects of the magnetic field in the kinetic and 
Rashba Hamiltonians via their dependence on quasiparti- 
cle momentum. The electron and hole degrees of freedom 
are represented, as usual, with Pauli matrices T x y z acting 
in the so-called Nambu space. 

The kinetic energy is split into non-magnetic and mag- 



netic terms as Hkin = 



■^*2> where 



H 



(0) 
kin 



2m* 



(1) 



H 



(i) 

/, / u 



2m* 



£2±, 



Ml 



J . ,Pz n P , Pz . P . 2 i 

P 2 2 - o P 5 

cos 6 t z 



02 02 



sin cos ( 



(2) 



We employ the cylindrical coordinate system indicated in 
Fig. [T] and denote the chemical potential by p, and the 
magnetic length along a general direction a = x,y 7 z by 
£ a = ^Jhc/eB a . Notice that, trivially, when a certain field 
component vanishes the corresponding magnetic length di- 
verges, giving a vanishing contribution to n^ n - 

The Rashba Hamiltonian originates in an effective elec- 
tric field £ as Hr oc a ■ {p x £ ) . We initially considered 
a radial field £ = £u p but, since this does not lead to the 
appearance of any Majorana states, we finally assumed a 
uniform direction £ = £u x . Similarly to the kinetic term, 
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The Zeeman contribution reads 
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where n gives the direction of the magnetic field and pa- 
rameter Ab is related to the gyromagnetic factor g by 
Ab = gfJ-BB/2. The final Hamiltonian contribution is the 
superconductivity (pairing) one 
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In summary, the Hamiltonian depends on the field ori- 
entation n, as well as on parameters Ao, Ag and a rep- 
resenting the superconducting, Zeeman and Rashba cou- 
pling strengths, respectively. In order to investigate realis- 
tic setups we assume parameter values from Ref . [24] . The 
effective mass, Rashba spin-orbit coupling, and radius of 
InSb nanowire are m* = 0.015m e (m e is the electron bare 
mass), a ~ 0.2 eVA and p w 50 nm. A convenient choice 
of energy and length units is then 



E V = 
L v = 



n 2 



= 2.03 mcV 



p = 50 nm 
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Fig. 2: Eigenvalues vs As when orbital effects are omitted. 
Panels (a), (b) and (c) correspond to applied magnetic fields 
pointing along x, y and z, respectively, x and z directions are 
almost identical. Parameters: A = 0.12Eu, a = 0.2EuLu, 
H — OEu, p = Lu, and L z = 30Lu. Number of basis: N n = 5f 
and N m = 31. Only the twelve eigenvalues closer to zero energy 
have been displayed. 

We also consider a superconducting parameter Ao « 0.25 
mcV and a Zceman energy that for g « 50 is related to 
the field B by A B = 0.723BT- 1 ^ [/ - 

Finite cylinders. — We have performed numeri- 
cal diagonalizations of the Hamiltonian for a cylinder 
of length L using as a basis states \nms a s T ) , where 
n = 1, 2, . . . represent square well eigenstates (z motion), 
m = 0, ±1, ±2, . . . represent L z angular momentum eigen- 
states {<f> motion) while s a = ± and s T = ± are the usual 
spin and isospin two component eigenstates. 

Figures [5] and [3] show the energy spectrum of a long 
cylindrical nanowire (L = ZOLjj) when the magnetic field 
is applied along the three Cartesian axis. In Fig. [5] orbital 
effects have been neglected (£ x . y . z = oo) while in Fig. [3] 
the full Hamiltonian is considered. For the y direction the 
spectrum in Fig. [5p is dense, with accidental crossings at 
zero energy but no robust Majorana modes. Field orienta- 
tions along x and z (Figs.[2ji and [2};) are basically equiva- 
lent and they show the emergence of a zero energy (Majo- 
rana) mode beyond a critical field A^' = 0.15Eu. Right 
after the critical field a clear energy gap separates this 
mode from other excitations. However, as Ab increases 
an enhanced oscillating behavior manifests the degrada- 
tion of the Majorana mode. This is a finite size effect and 
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Fig. 3: Same as Fig. f2]with orbital effects. 

agrees with our previous study of planar geometries |29j 
as well as with Ref. [30) . Importantly, recent experimental 
evidences have confirmed this remarkable dependence of 
the Majorana pair splitting by measurements of the tun- 
neling differential conductance [28]. When Ab exceeds 
0.5Ejj two more modes emerge in Fig. [2^ and Et; i.e., 
a total of three Majorana modes are oscillating around 
zero energy. Labeling each topological phase of the cylin- 
der by the number of different Majorana modes the ex- 
pected sequence when increasing the magnetic field is then 
— 1 — 3 — ... . Each m-subset contributes an indepen- 
dent Majorana after some critical field, while degeneracy 
of m and — m causes modes to emerge in pairs, except for 
m = 0. We focus next on the changes to this scenario 
induced by orbital effects. 

Figure [3] shows the influence of orbital effects on the 
nanowire energy spectrum. Now x and z orientations are 
no longer equivalent and a large anisotropy in the gap- 
closing field is obtained. This critical field is strongly re- 
duced along the sequence x — > y — > z. Along x (Fig. 
[3K) a state of unclear Majorana character emerges beyond 
Ab ~ 0.5-E[/, with an oscillating energy and nearby ex- 
cited modes. For the z direction (Fig. \3jp) a clear Ma- 
jorana mode emerges at a much lower magnetic field, 
Ab ~ Q.\hEij. Interestingly, when increasing the field 
along z a second Majorana mode seems to be present 
for 0.18-E[/ < A B < 0.27E L r. This behavior is in sharp 
contrast with planar geometries, where Rashba mixing 
precludes the coexistence of more than one Majorana 
mode. In cylinders, Rashba mixing is effectively weaker, 
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Fig. 4: Density distributions corresponding to the two Majo- 
rana modes of Fig. [3J; with A_b = 0.2Eu- An arbitrary gray 
colour scale is used, with dark colour indicating higher density. 

as its strength changes sign with the azimuthal angle [cf. 
Eqs. <j3j> and Q]. Notice also that the spin-orbit length, 
Lso = h 2 /ma « 250 nm, is Lso >> Pi also indicating a 
regime of weak spin-orbit where multiple Majoranas may 
coexist [31]. Increasing further the magnetic field in Fig. 
[3H, oscillations around zero energy become large and the 
Majorana character of the states is blurred. The sequence 
of topological phases as the field is increased suggested in 
Fig. [3J; is 0-1-2-1-2-3-2. This is in clear difference 
with the findings in absence of orbital effects, Fig. [2J;. 

In transport measurements, the presence of multiple 
Majorana modes affects the zero-bias conductance Go- 
More precisely, Go = MQ, where M is the number of 
Majorana modes and Q is a quantization unit depending 
on the nature of the contacts (Q = 2e 2 /h for a normal- 
superconductor junction). Although, currently, experi- 
ments are still far from this limit of high conductances 
|32j , the observation of a zero bias conductance following 
the above sequence with magnetic field would be a clear 
proof of the cylinder topological phases. 

As an illustrative case, Fig. 2] shows the densities of the 
two Majorana modes of Fig. [3): when = 0.2Ejj- As 
expected, the densities are localized close to the cylinder 
edges and decay towards the central part. Remarkably, 
however, the density maxima of the two modes take com- 
plementary azimuthal positions, indicating an effective re- 
pulsion between Majorana modes. We have also obtained 
that when the mode energies deviate from zero due to the 
finite size effect, for large Ag's, the densities extend more 
towards the centre, eventually occupying all of the cylin- 
der length. 

Semi-infinite cylinders. — The above results show 
that finite size effects are relevant for an aspect ratio 
L/p = 30. We have checked that for L/p = 60 (re- 
sults not shown) the finite size oscillations are reduced 



but the behavior is qualitatively the same at high fields. 
Therefore, in order to ascertain the sequence of topological 
phases mentioned above, we have directly addressed the 
semi-infinite cylinder using the complex band structure 
approach of Ref . [33J . In the semi-infinite system there is 
a single edge, the other one lying at infinite distance, and 
finite-size artifacts are absent by definition. 

For a complex wave number k the eigenstate of the 
infinite-cylinder Hamiltonian reads 

^,^V„Vr) = e ikz J2^BMXsAv,)XsAVr) , (8) 

where, generically, rj a _ T =f, \. refers to a spin-isospin dis- 
crete variable and s a , T = ± is the corresponding quantum 
number. We have determined the allowed wave numbers 
k and state amplitudes V'saSr f° r a given energy using the 
approach of Ref. [33] . A Majorana mode is signaled by a 
zero eigenvalue in matrix 

M k « = Y,Jd<t> 4% (0) (0) • (9) 

In practice the set of complex wave numbers is truncated 
to a finite size using a cut-off in momentum and the emer- 
gence of a zero eigenvalue of matrix M is actually seen 
as a steadily decreasing eigenvalue with increasing value 
of the cut-off. The approach of Ref. [33] provides also a 
method to obtain the critical fields corresponding to phase 
transitions. They are seen as zeros in function JF(Ab), the 
derivative discontinuity at an arbitrary matching point for 
k = 0. 

Figure [5] confirms the scenario inferred above from the 
finite cylinder calculations. Remarkably, the sequence of 
nodes in T (Fig. [5h) has a clear correspondence with the 
changes in the finite cylinder spectrum of Fig. [3J:. In 
the semi-infinite cylinder the critical fields for topologi- 
cal transitions can be determined accurately, while in the 
finite cylinder transitions are blurred, specially at high 
fields. Lower panels in Fig. [5] contain the evolution with 
cut-off in wave numbers of the eigenvalues of matrix Ai. 
In the first region (A B < 0.15Eu) the lowest eigenvalue 
converges to ~ 1, indicating that no edge mode is possi- 
ble. In the second and third regions we see that one and 
two eigenvalues, respectively, keep decreasing as the cut- 
off is increased. This is a clear evidence that second and 
third regions contain one and two Majorana modes, re- 
spectively. Similar plots (not shown) are obtained for the 
successive regions, containing the number of Majoranas 
indicated in Fig. [5^. Topological invariance is manifested 
as independence of the number of Majorana modes with 
the precise value of A#, as long as this remains within a 
given interval between nodes of Fig. EK- 

The spatial densities of the two Majorana modes for 
Ab = 0-2Ejj are shown in Fig. [6] When compared with 
the finite cylinder, densities are in good agreement for 
cases when the finite-cylinder mode is very close to zero 
energy. If not, large discrepancies in the bulk of the 
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Fig. 5: Upper panel: function .F whose nodes signal the topo- 
logical transitions. The number of Majorana edge states are 
indicated for the low field regions. Lower panels: Evolution 
of matrix A4 eigenvalues with the cut-off in wave numbers (in 
absolute value). Each eigenvalue tending to zero in the limit 
of large |fc|'s corresponds to a Majorana mode. Except for the 
infinite cylinder length all other parameters are the same of 
Fig. [2b. 



cylinder are found. An obvious difference, of course, is 
that in the finite cylinder the density is repeated on the 
two edges. A similar avoided mode overlap is found in 
both cases, and there is only a minor angular shift in 
the relative position of the maxima. We finally men- 
tion that the wave function of Majorana modes is quasi- 
real, i.e., its imaginary part almost vanishes. The wave 
function is more real the more fe's are included in its 
expansion, and fulfills &(z, <f>, t> f) = — ^{z, 4>i\-iV) an d 
^{z, <fi, f, i) = ^(z, 4>, i, t)- Both things are to be expected 
for proper Majorana solutions. 

Conclusions. — In closing, we have investigated the 
occurrence of Majorana zero-energy modes at the edges 
of cylindrical nanowires using both finite and semi-infinite 
models. We have considered a Rashba field along x and 
a magnetic field in arbitrary direction. Under this situa- 
tion we report the strong influence of the orbital effects 
on the appearance of zero energy states. First, the critical 
fields corresponding to gap closings strongly depend on 
the magnetic field orientation. Second, in contrast with 
planar geometries we observe phases where several Majo- 
rana modes coexist. Orbital effects define the sequence 
of phases with increasing field. The observation of the se- 
quence 0—1 — 2 — 1 — 2 — 3... in the zero bias conductance 
as a function of magnetic field (in units of the conductance 
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Fig. 6: Density distributions corresponding to the two Majo- 
rana modes of a semi-infinite cylinder with As = 0.2Eu- All 
other parameters are the same of Fig. [2b. 



quantum) would be a clear manifestation of the cylinder 
topological phases. As in planar geometries, at high fields 
an oscillating behavior signals the degradation of the Ma- 
jorana character due to finite size effects. 
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